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Simulation  Subsection 


GOAL 

For  Phase  I,  the  goal  of  this  project  was  to  establish  a  protocol  or  pseudo-automated  procedure  for  classical 
force  field  parameterization.  The  force  field  takes  the  form: 
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Parameters  required  include  kb,  r0,  ke,  90,  kz,  n,  8,  k,|(,  \j/0,  Sq,  rminjij,  q,,  qj.  A  combination  of  ab  initio  electronic 
structure  calculations  (using  GAMESS)  and  molecular  mechanics  calculations  was  used  to  optimize  parameters 
for  a  user-specified  molecule. 


PHASE  I  PROGRESS/RESULTS 
Test  molecules: 

(A)  1 -Ethyl- 1-methylphospholinium  Tricyanomethanide 

(some  parameters  for  1 -ethyl- 1-methylphospholinium  cation  taken  from  the  work  of  Canongia  Lopes 
and  Padua)  1 


(B)  1  -Ethyl-3 -methylimidazolium  Butanesulfonate 

(some  parameters  for  1 -ethyl-3 -methylimidazolium  cation  taken  from  the  work  of  Cadena  et  al).2 


I.  Intramolecular  terms  (bond,  angles,  dihedrals,  impropers) 


I.a.  Bonds  and  Angles 

The  GAMESS3  source  code  has  been  modified  to  punch  optimized  geometry  parameters  (r0,  0o)  and  the  force 
constant  matrix  (containing  the  force  constants  kb  and  ke)  for  a  given  molecule.  The  force  constant  matrix  can 
be  punched  in  internal  coordinates  (with  appropriate  units  of  Hartree/Bohr~  for  bonds  and  Hartree/rad '  for 
angles)  according  to  how  the  molecule  was  designated  in  the  Z-matrix  (not  unique).  A  protocol  is  used  to  write 
‘intrinsic  frequencies’  corresponding  to  localized  vibrational  modes, 4  allowing  one  to  extract  intrinsic  force 
constants  for  bond  stretching  only.  These  can  be  read  directly  into  a  fonnatted  (i.e.,  CHARMM)  force  field. 
For  the  test  molecules  1 -Ethyl- 1-methylphospholinium  Tricyanomethanide  and  1 -Ethyl-3 -methylimidazolium 
Butanesulfonate,  the  geometries  were  optimized  at  the  B3LYP/TZVP  level5,  and  the  force  constants  were  used 
at  this  level,  with  appropriate  scaling  based  on  small  molecule  comparisons. 6 

Due  to  strong  off-diagonal  couplings  in  angle  bending  modes,  the  intrinsic  force  constants  for  angle  bends  will 
be  more  difficult  to  acquire,  as  their  magnitudes  in  the  force  constant  matrix  are  generally  too  high.  This  is 
evidenced  in  the  example  of  the  heterocyclic  cation  1 -butyl-3 -methylimidazolium,  for  which  the  diagonal  force 
constant  element  for  the  N-C-N  ring  angle  is  punched  to  be  ~3  Hartree/rad“,  with  several  off-diagonal  couplings 
of  the  same  magnitude.  This  may  be  rationalized  due  to  restrictive  ring  constraints.  This  issue  is  also  relevant 
to  other  systems  where  dependent  angle  bending  modes  exist.  A  solution  is  required  for  this  problem,  such  as 
making  the  force  constant  matrix  as  diagonal  as  possible,  thus  generating  the  most  localized  force  constants 
possible.  In  addition,  judicious  choice  of  internal  coordinates  is  required,  as  different  initial  Z-matrices  or 
poorly  chosen  coordinates  may  produce  different  or  erroneous  results.  For  most  systems,  the  diagonal  angle 
force  constants  may  possibly  be  used,  with  the  qualification  that  they  are  probably  an  upper-bound  to  the  true 
values.  The  quality  of  the  force  constants  used  must  necessarily  go  hand-in-hand  with  the  desired  result.  If  the 
calculated  property  is  only  minimally  dependent  on  the  intramolecular  parameters  (i.e.,  static  properties),  then 


the  amount  of  effort  imparted  to  achieve  the  desired  outcome  must  be  evaluated.  Due  to  these  challenges,  the 
angle  force  constants  were  derived  by  performing  structural  perturbations  on  smaller  related  molecules  where 
the  angles  were  incremented  by  -2  to  2°,  followed  by  single  point  energy  calculations  at  the  B3LYP/TZVP 
level.  The  force  constants  were  derived  assuming  a  harmonic  approximation  and  averaged. 

GAMESS-specific  required  flags  and  brief  description: 

Sstatpt  hssend  =  .t.  (compute  hessian  at  stationary  point) 

Sforce  vibanl  =  .t.  (perform  vibrational  analysis) 

decomp  =  .t.  (activates  internal  coordinate  frequency  analysis) 

prtifc  =  .t.  (print  internal  force  constants) 

puncff  =  .t.  (print  force  field  information  at  end  of  log  file) 

(Note:  The  puncff  flag  was  specifically  added  for  this  project). 

Lb.  Dihedrals 

Dihedral  analysis  and  parameterization  requires  a  combination  of  ab  initio  calculations  and  molecular 
mechanics  calculations.  A  perl  script  has  been  written  ( torsion.pl)  to  generate  a  series  of  GAMESS  input  files 
for  dihedral  scan  analysis.  The  optimized  geometry  for  a  given  molecule  (performed  previously  in  section  I.a.) 
is  placed  into  a  GAMESS  input  file  (i.e.,  example. inp).  The  geometry  must  be  specified  in  Z-matrix  form. 
Additionally,  the  user  must  set  a  flag  corresponding  to  the  particular  dihedral  that  will  be  scanned  and  held 
fixed  during  the  constrained  optimization. 

GAMESS-specific  required  flag  needed  for  the  perl  script  (and  description): 

Sstatpt  ifreeze(  1)  =  x  (x  is  numerical  reference  to  the  particular  internal 

dihedral  coordinate  to  be  scanned  as  defined  in  Z-matrix) 

The  perl  script  is  executed  as  follows  for  the  example. inp  case  file: 

./torsion.pl  example. inp  y  example. d#### 

The  variable  y  indicates  the  exact  line  of  the  Z-matrix  containing  x  (noting  that  all  lines  containing  flags  are 
ignored  by  the  script).  The  perl  script  creates  a  series  of  input  files  (containing  the  flags  listed  in  example. inp) 
in  which  the  particular  dihedral  (x)  to  be  scanned  in  the  molecule  has  been  rotated  in  +10°  increments.  The  files 
created  are  of  the  type  as  shown  below,  where  ####  are  the  specific  atoms  labels  in  the  dihedral: 

example. d####+10. inp,  example. d####+20. inp,  ...  ,  example. d####+360. inp 

If  a  series  of  unique  dihedral  angles  are  present  in  the  molecule,  individual  example. inp  files  must  be  created 
with  the  unique  ifreeze(l)  numerical  references  specified,  and  the  perl  script  is  run  for  each  of  these  to  create 
the  scan  files.  A  system-specific  batch  submission  script  may  be  written  to  run  these  GAMESS  input  files. 

It  is  general  practice  in  force  field  development  to  run  constrained  optimization  scans  at  a  low  to  moderate  level 
of  theory.  Following  the  scan,  the  optimized  geometries  at  each  scan  point  must  then  be  extracted  and 
reinserted  into  new  GAMESS  input  decks  in  which  the  theory  level  has  been  modified  to  obtain  proper 
energetics  as  a  function  of  dihedral  angle.  For  the  anions  of  test  molecules  of  (A)  and  (B),  appropriate  dihedral 
scans  were  performed  at  the  B3LYP/TZVP  level  of  theory.  It  was  found  that  these  profiles  closely  matched 
those  calculated  at  the  MP2(FC)/cc-pVTZ(-f)  level,  where  (-/)  means  that  f-type  functions  have  been  removed 
from  the  basis  set.  7 

Following  the  ab  initio  analysis,  a  molecular  mechanics  (MM)  analysis  is  also  required.  Here,  the 
B3LYP/TZVP  optimized  geometries  are  also  extracted  and  inserted  into  .pdb  files  (protein  databank  files)  or 
other  coordinate  files.  The  .pdb  files  are  commonly  used  in  various  MM  packages  such  as  CHARMM8  or 


NAMD9.  Connectivity  must  also  be  specified,  and  tools  exist  to  create  proper  connectivity  files  (as  well  as  to 
create  .pdb  files).  From  here,  relaxed  MM  energies  may  be  obtained  using  a  force  field  in  which  the  particular 
dihedral  parameters  begin  studied  are  set  to  zero  and  a  constrained  MM  minimization  is  performed.  It  is  critical 
that  all  other  force  field  parameters  (for  bonds,  angles,  charges,  and  Lennard-Jones  interactions)  be  pre¬ 
specified.  The  MM  energies  are  then  subtracted  from  the  ab  initio  energies,  and  the  resultant  curve  is  then  fit  to 
a  cosine  series  using  non-linear  regression  tools.  The  number  of  terms  in  the  cosine  series  is  varied  to  obtain 
the  best  possible  fit  the  data. 

Challenges  that  still  remain  for  obtaining  dihedral  force  field  parameters  mainly  involve  automating  the  step¬ 
wise  procedure  of  ab  initio  scan,  MM  scan,  subtraction,  and  non-linear  regression  as  described  above. 
However,  all  of  the  tools  are  in  place  to  reach  the  goal  and  proof-of-concept  may  be  demonstrated  with  a  well- 
defined  procedure,  and  automation  of  this  step  may  be  implemented  at  a  later  date.  A  significant  challenge 
remains  for  the  treatment  of  ring  systems  (particularly  ionic  liquids),  since  scanning  of  the  ring  dihedrals  may 
produce  unphysical  results.  In  this  case,  identifying  the  symmetry  of  the  problem  and  tailoring  the 
parameterization  accordingly  may  be  critical. 

For  Phase  I,  the  above  described  procedure  was  performed  manually  to  obtain  dihedral  parameters  for  the 
anions  of  (A)  and  (B).  Dihedral  parameters  for  the  cation  of  (B)  were  taken  directly  from  the  work  of  Cadena  et 
al.2  Dihedral  parameters  for  the  cation  of  (A)  were  derived  from  the  work  of  Canongia  Lopes  and  Padua  1  for 
linear  phosphonium  cations  and  applied  to  the  cyclic  phospholinium  system. 

II.  Intermolecular  terms  (Point  charges,  Lennard-Jones  parameters) 


II. a.  Point  Charges 

As  atomic  centered  charges  are  not  quantum-mechanical  observables,  many  procedures  exist  to  assign  these 
charges  from  ab  initio  electronic  structure  calculations.  Force  field  charges  are  primarily  assigned  using  the 
CHELPG  procedure,  10  in  which  the  electrostatic  potential,  calculated  at  a  particular  level  of  theory,  is  projected 
onto  the  atomic  sites,  such  that  the  gas-phase  dipole  moment  of  the  molecule  is  conserved.  Charges  for  the 
cations  and  anions  of  (A)  and  (B)  were  calculated  at  the  MP2/cc-pVTZ(-f)//B3LYP/TZVP  level  of  theory.  11 
GAMESS  is  equipped  to  routinely  calculate  CHELPG  charges. 


GAMESS-specific  required  flags  and  brief  description: 

Selpot  iepot  =  1  (initialize  electrostatic  potential  calculation) 

where  =  pdc  (type  of  calculation) 

$pdc  ptsel  =  chelpg  (calculate  chelpg  charges) 

constr  =  dipole  (constraint  to  reproduce  the  dipole  moment) 


(Note:  The  puncff  flag  described  above  will  punch  the  chelpg  charges  at  end  of  file). 


II. b.  Lennard-Jones  Parameters 

Lennard-Jones  (LJ)  parameter  assignment  and  optimization  poses  the  most  significant  challenge  to  an 
automated  force  field  generation  protocol.  Moderately  complex  procedures  exist  to  obtain  LJ  parameters  from 
ab  initio  electronic  structure  calculations.  For  example,  an  atom  with  well-known  LJ  parameters  (i.e.,  Ar)  is 
moved  along  the  surface  of  a  molecule,  and  at  particular  locations,  the  interaction  energies  are  calculated  (high 
level  calculations  are  required).  Assuming  a  given  set  of  mixing  rules  and  a  defined  mathematical  form  for  the 
interaction  potential,  the  LJ  parameters  may  be  obtained  using  appropriate  fitting  procedures.  Automation  of 
the  LJ  parameterization  using  a  similar  protocol  is  an  appropriate  goal  for  Phase  II/Phase  III  of  the  project. 
Currently,  the  protocol  involves  extracting  LJ  parameters  from  predefined  force  fields  (i.e.,  CHARMM,  Amber) 
or  literature  sources  for  the  atoms  in  similar  chemical  environments  as  exist  if  the  user-defined  molecule.  For 
example,  an  ionic  liquid  containing  a  ether  group  in  the  alkyl  chain  potion  of  the  molecule  would  require  LJ 
parameters  for  oxygen  in  a  similar  environment  as  described  in  another  force  field.  Due  to  the  vast  amount  of 
force  field  data  that  exists  in  the  scientific  community,  ample  opportunity  is  present  to  tabulate  much  of  the  data 


and  obtain  the  best  possible  set  of  LJ  parameters  for  a  diversity  of  atom  types  to  be  used  in  Phase  I  of  the 
automation  procedure.  We  are  currently  developing  a  force  field  database  from  which  appropriate  LJ 
parameters  can  be  found  from  compiled  literature  sources.  LJ  parameters  for  the  anions  were  taken  from  the 
CHARMM  force  field8  for  similar  atom  types,  and  the  LJ  parameters  for  the  cation  were  taken  from  literature 
sources.  1 2  A  beta  version  of  the  database  exists  at:  http://www.nd.edu/~hnmr 

GAMESS  possess  the  capability  to  calculate  the  Lennard- Jones  parameters  using  an  entirely  first  principles 
approach,  thus  allowing  for  the  development  of  internally  consistent  values  for  force  field  work  and  eliminating 
the  need  to  pick  and  choose  literature  values.  Here,  an  effective  fragment  potential  (EFP)  method  is  used  to 
develop  interaction  potentials  based  on  localized  orbitals.  Thus,  tools  are  available  for  the  parameterization 
using  electronic  structure  theory,  and  require  projection  of  the  localized  orbital  values  onto  the  atoms  to  achieve 
the  desired  parameters.  This  approach  is  considered  to  be  a  Phase  11/Phase  III  objective. 

It  is  noted  that  force  field  parameters  may  be  refined  in  particular  cases  based  on  comparisons  to  available 
experimental  data.  For  example,  comparison  of  simulated  liquid  densities  or  heats  of  vaporization  to 
experimental  values  may  warrant  refinement  of  the  rm;njij  or  eij  values.  These  properties  are  to  a  large  extent 
decoupled  and  therefore  the  LJ  parameters  may  be  refined  fairly  independently  of  one  another.  It  is  also  noted 
that  this  refinement  requires  the  user  to  evaluate  the  quality  of  the  simulation  data,  and  involves  a  greater 
amount  of  user  intervention  than  would  be  required  for  the  general  force  field  development  procedure  described 
herein. 

III.  Other  Considerations 

The  above  outlined  strategic  approach  to  force  field  parameterization  will  be  immensely  beneficial  to  the  ionic 
liquid  community  in  the  development  of  new  ionic  liquids  for  specific  tasks  and  the  evaluation  of  their 
properties.  The  tremendous  number  of  possible  ionic  liquids  that  are  within  reach  makes  it  critical  that  a  reliable 
approach  is  available  to  guide  experimental  efforts  and  design.  In  addition,  the  automation  of  force  field 
development  will  also  be  beneficial  to  other  fields,  as  the  approach  in  not  restricted  to  ionic  liquids.  The  ability 
to  have  reliable  force  field  at  one’s  fingertips  for  the  study  of  materials  that  exhibit  non-linear  optical  properties 
or  for  materials  that  contain  transition  metals,  lanthanides,  and  actinides  will  also  have  a  positive  impact  in  the 
simulation  community. 


PHASE  I  PROTOCOL  OVERIEW:  CONCLUSIONS 


1 .  DEFINE  AND  DRAW  MOLECULE  TO  BE  PARAMETERIZED 

a.  Use  wxMacMolPlot  (version  7.0b3) 

b.  Define  a  Z-matrix  for  the  molecule.  Use  a  Z-matrix  that  captures  the  largest  number  of  required 
parameters 

c.  Assign  individual  atom  types  consistent  with  force  field  development 

2.  OPTIMIZE  EQUILIBRIUM  GEOMETRY  AND  CALCULATE  VIBRATIONAL  FREQUENCIES 
OF  MOLECULE  USING  GAMESS 

a.  Use  consistent  level  of  theory,  as  common  in  force  field  development 

b.  Screen  for  imaginary  modes;  if  present,  modify  geometry  and  try  again 

c.  Use  appropriate  GAMESS  flags  to  punch  force  field  infonnation 

3.  EXTRACT  BOND  DISTANCES  AND  ANGLES  AND  CORRESPONDING  FORCE  CONSTANTS 
FROM  GAMESS  OUTPUT 

a.  Use  intrinsic  force  constants  for  bonds 

b.  Use  diagonal  elements  of  internal  force  constant  matrix  for  angles 

c.  Check  for  presence  of  ring,  and  assign  force  field  parameters  accordingly 

d.  Repeat  steps  1-3  for  internal  coordinates  that  were  missing  from  original  Z-matrix,  since  not  all 
internals  are  necessarily  uniquely  defined 

e.  Export  values  to  force  field  parameter  file  (.PRM)  in  appropriate  fonnat  (i.e.,  CHARMM) 

4.  EXTRACT  CHELPG  ATOMIC  CHARGES  FROM  GAMESS  OUTPUT 

a.  Compare  CHelpG  dipole  with  ab  initio  dipole  to  ensure  agreement 

b.  Export  values  to  force  field  connectivity  file  (.PSF)  in  appropriate  format  (i.e.,  CHARMM) 

5 .  ASSIGN  LJ  PARAMETERS  FROM  TABULATED  DATA 

a.  Read  atoms  types  defined  in  step  lc  and  identify  corresponding  LJ  parameters 

b.  Export  values  to  force  field  parameter  file  (.PRM)  in  appropriate  fonnat  (i.e.,  CHARMM) 

c.  Use  GAMESS  EFP  method  when  available 

6.  PERFORM  AB  INTIO  DIHEDRAL  OPTIMIZATION  SCANS  AND  ENERGETICS  SCANS 

a.  Create  GAMESS  input  file  with  appropriate  flags  and  Z-matrix  corresponding  to  optimized 
geometry  obtained  in  2a. 

b.  Identify  dihedral  angle  to  scan  (internal  coordinate  number  and  Z-matrix  line  number) 

c.  Run  torsion.pl  script  on  this  input  file  to  create  series  of  GAMESS  scan  files  (36  total) 

d.  Use  batch  submission  script  to  run  each  GAMESS  constrained  optimization  calculation  where 
assigned  dihedral  coordinate  is  kept  fixed 

e.  Extract  optimized  coordinates  from  each  of  the  36  scan  files  and  export  into  36  new  GAMESS 
input  files  for  high  level  B3LYP/TZVP  or  MP2/cc-pVTZ(-f)  energy  calculations 

f.  Use  batch  submission  script  to  run  each  GAMESS  single  point  energy  calculation  where 
geometry  is  kept  fixed 

g.  Plot  energy  vs.  dihedral  angle  to  observe  periodicity  and  check  for  problems 

7.  PERFORM  MOLECULAR  MECHANICS  DIHEDRAL  RELAXATION  SCANS  (EXAMPLE  W/ 
NAMD) 

a.  Extract  optimized  coordinates  from  each  of  the  36  scan  files  (6d)  and  export  into  36  new  .PDB 
files  for  MM  energy  calculations 

b.  Define  molecular  connectivity  in  .PSF  file 

c.  Define  MM  simulation  details  in  .CONF  file,  including  temperature  (0  K),  constraints  (dihedral), 
and  minimization  steps 


d.  Confirm  that  force  field  file  (.PRM)  is  complete,  with  dihedral  parameters  set  to  zero  for  dihedral 
being  studied 

e.  Use  batch  submission  script  to  run  each  NAMD  minimization  calculation,  where  all  other 
degrees  of  freedom  allowed  to  relax 

f.  Plot  energy  vs.  dihedral  angle  to  observe  periodicity  and  check  for  problems 

8.  PERFORM  SUBTRATION  AND  COSINE  SERIES  FIT 

a.  Subtract  curve  7f  from  curve  6g  to  obtain  energy  vs.  dihedral  angle  curve  to  be  fit  with  cosine 
series 

b.  Use  nonlinear  regression  to  obtain  best  fit  cosine  series  parameters,  where  the  number  of  terms  in 
the  series  is  altered  to  obtain  ideal  fit 

c.  Export  values  to  force  field  parameter  file  (.PRM)  in  appropriate  fonnat  (i.e.,  CHARMM) 

d.  Repeat  steps  6  to  8  for  all  dihedrals  to  be  fit,  keeping  in  mind  that  the  ideal  parameterization 
procedure  may  be  a  ‘building-up’  procedure  (Padua,  J.  Phys.  Chem.  A  2002,  16,  10116),  starting 
with  smaller  molecules  of  similar  chemical  composition  OR  may  require  appropriate  dihedral 
constraints/ restraints . 

9.  RUN  SIMULATIONS  AND  REFINE  FORCE  FIELD  IF  NECESSARY 

a.  Run  molecular  simulation  (Molecular  Dynamics  or  Monte  Carlo)  using  newly  parameterized 
force  field. 

b.  Compare  several  calculated  properties  (i.e.,  liquid  density,  heat  of  vaporization)  with  available 
experimental  data  and  evaluate  quality  of  force  field 

c.  If  necessary,  refine  the  LJ  parameters  and  repeat  steps  5-8. 


FORCE  FIELDS  DEVELOPED  FOR  TEST  ILs  (A)  and  (B) 


(A)  l-Ethyl-l-methylphospholinium  Tricyanomethanide 
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(B)  l-Ethyl-3-methylimidazolium  Butanesulfonate 
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CALCULATED  PROPERTIES 


Results 


(A)  1  -Ethyl-1  -methyl-phospholinium  Tricyanomethanide 


Liquid  density  (calculated) 

p  (g/cm3):  1.062 


Self-diffusion  coefficient  (calculated) 

Ds  (cation,  m2/s)  :  1.08  x  10'L’ 

Ds  (anion,  m2/s)  :  1.72  x  10‘13 


(B)  l-Ethyl-3-methylimidazolium  Butanesulfonate 

Liquid  density  (calculated) 

p  (g/cm3):  1.133 

Self-diffusion  coefficient  (calculated) 

Ds  (cation,  m2/s)  :  2.72  x  10'12 
Ds  (anion,  m2/s)  :  1.45  x  10‘12 

Experimental  Subsection 

A.  Project  Justification  and  Synthetic  Target  Selection: 

Hydrazinium  Salts  as  Potential  Energetic  Material  Candidate: 

Hydrazinium-based  ionic  liquids  (ILs)  have  been  selected  as  potential  candidates  for  use  in  propulsion 
applications,  as  it  has  been  demonstrated  previously  that  such  compounds  can  be  successfully  utilized  as 
propellants.  '  Additionally,  as  ionic  liquids,  these  compounds  are  found  to  bear  other  favorable  properties 
that  are  less  frequently  observed  for  materials  of  strictly  molecular  (non-ionic)  nature:  (i)  relative  insensitivity 
to  frictional  stress  and  electrostatics,  (ii)  negligible  volatility  due  to  very  low  vapor  pressures,  and  (iii) 
relatively  good  thennal  stability.  Additionally,  due  to  the  ionic  liquid  nature  of  these  materials  they  have 
advantage  of  broad  liquidus  ranges  and  the  potential  for  dual-functionalization  of  either  component  anion  or 
cation. 

2  -Hydroxyethyl  -Hydrazinium  Salts:  Justification  for  the  choice  of  the  specific  cation  for  investigation: 

Salts  based  on  the  parent  structure  of  the  neutral  2 -hydroxy ethyl-hydrazine  (HEH)  have  been  previously 
reported  by  Brand  and  Drake,  where  it  was  anticipated  that  the  addition  of  the  2 -hydroxyethyl  side-chain  would 
result  in  reduced  volatility  in  comparison  to  hydrazine. 14  Such  salts  have  been  chosen  as  the  subjects  for  a 
computational  investigation,  and  the  objective  for  the  team  here  at  The  University  of  Alabama  is  to  serve  as  the 
center  for  the  synthesis  of  chosen  targets  and  to  provide  all  needed  data  that  will  support  the  computationally 
modeled  systems. 

2-Hydroxyethyl-Hydrazinium  Salts:  Targeted  1:1  and  1:2  salts 

The  mono-  and  diprotic  salts  of  hydrazinium  have  been  selected  as  specific  mono-  and  dianionic  targets, 
respectively,  and  they  are  obtained  by  protonation  of  HEH  with  a  broad  array  of  common  Bronsted  acids 
(Figure  1): 


hn-nh2 


HA  ini  :1  or  1:2  equiv. 
r.t.  24  h 


HO 


© 

h2n-nh2 


OR 


©  © 

h2n-nh3 


a  =  no3',  cr,  bf4',  pf6",  n(so2cf3)' 

Figure  1.  Mono-  and  diprotonation  of  2-hydroxyethyl-hydrazine  by  addition  of  common  Bronsted  acids 


♦  The  synthetic  targets  generated  will  permit  comparison  of  change  in  anion  identity  as  well  as 
comparison  of  mono-  and  dicationic  salts. 

♦  Perchlorate  (GOT)  was  originally  indicated  as  a  potential  target,  but  it  has  been  removed  from  the 
current  set  due  to  potential  instability  issues  in  handling  the  material. 


B.  Experimental  Section: 

Work  Plan  and  Status: 


♦  Synthesize  2-hydroxyethyl-hydrazinium  salts  in  1:1  and  1:2  ratios  (catiomanion)  as  ~  1.5  gram  sample: 

o  Nitrate  (NOT)  [DONE] 
o  Chloride  (Cf)  [DONE] 
o  Tetrafluorob orate  (BF4‘)  [DONE] 
o  Hexafluorophosphate  (PF6‘)  [DONE] 
o  Bistrifylamide  (NT1V)  [DONE] 
o  Perchlorate  (GOT)  [Target  postponed,  see  last  section] 

♦  Characterize  the  physical  properties  of  the  samples: 

o  1H-NMR  [DONE] 

o  Karl-Fischer  /  FT-IR  (determination  of  purity) 
o  TGA:  (Tonset,  T5o/odec,  T(|cc)  [DONE] 
o  DSC  (Tm,  Tg,  Tcryst,  etc.)  [DONE] 
o  Heat  capacity  (Cp) 

o  Thermal  Stability  (24  hour  @  75  °C  isothennal)  [DONE] 
o  XRD 


♦  Make  one  salt  in  5  gram  scale  (e.g.,  2-OHET-HY  [N03]  [DONE] 

♦  Order  2-OHEt-HY  from  Aldrich  [DONE] 

General  Synthetic  Approach  for  the  formation  of  2-hydroxyethyl-hydrazinium  (1:1)  salts: 

The  acids  were  used  as  concentrated  aqueous  solutions  in  reaction  with  2-hydroxyethyl-hydrazine 
(Fluka,  98%)  as  purchased.  To  a  20  mL  glass  vial,  the  2 -hydroxy ethyl  hydrazine  (10  mmol,  0.7610  g)  was 
added  followed  by  the  addition  of  acid  (10  mmol).  The  vials  were  capped  and  then  stirred  overnight,  followed 
by  ~12  hours  drying  in  air  at  50  °C  furnace. 

General  Synthetic  Approach  for  the  formation  of  2-hydroxyethyl-hydrazinium  (1:2)  salts: 

The  procedure  here  is  similar  to  that  described  above  for  the  1 : 1  salts,  except  that  the  amount  of  acid 
used  is  increased  to  20  mmol. 

General  Synthetic  Approach  for  the  formation  of  5  gram  scale-up  of  2-hydroxyethyl-hydrazinium  nitrate  (1:1) 
salt: 


The  synthesis  for  the  5  gram  scale  up  of  the  2-OHEt-HY[NC>3]  salt  is  similar  to  the  synthesis  of  other  1 : 1 
salts  described  above,  with  the  change  in  both  starting  materials  to  50  mmol  quantities  and  the  use  of  a  50  mL 
round-bottom  flask  as  the  reaction  vessel. 

X.  Results  and  Discussion: 

Part  A:  1 II-NMR 


Figure  Al:  'H-NMR  plot  for  2-hydroxyethyl  hydrazine  (neutral)  salt  precursor;  360  MHz  DMSO-</6 


♦  Highly  shielded  amine  and  hydroxy  proton  signal  (5  ~  3.865)  is  accurate  for  integration  of 
approximately  4H  =  (2HNh2  +  1HNh  +  IHoh). 

♦  -CH2OH  (8  =  3.430  ppm)  and  -CH2NH  (8  =  2.624  ppm)  appear  as  triplets  at  the  expected  integration 
and  chemical  shift  compared  with  literature  values. 


9.0  8.0  7.0  6.0  5.0  4.0  3.0 

ppm  (tl) 

Figure  A2:  'H-NMR  plot  of  2-hydroxyethyl  hydrazinium  mononitrate  (top)  and  dinitrate  (bottom)  salts;  DMSO-</6  360  MHz. 

♦  The  integration  and  down  Held  shift  of  amine/ammonium  and  hydroxy  protons  indicate  the  successful 
formation  of  the  mono-  and  dicationic  forms  of  the  nitrate  salts;  where  further  downfield  shifting  is 


observed  for  the  dicationic  salt  (2+  formal  charge,  and  more  deshielding  of  electron  density  from  the 
cationic  protons). 

♦  Small  1:1:1  triplet  at  5  ~  7. 1  ppm  is  indicative  of  the  Jn-h  spin-spin  coupling  that  is  typical  for  amine 
protons.  Further  Jhn-ch  coupling  is  not  observed  due  to  peak  broadening. 
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Figure  A3:  'H-NMR  plot  of  2-hydroxyethyl  hydrazinium  monochloride  (top)  and  dichloride  (bottom)  salts;  DMSO-</6  360 
MHz. 

♦  The  peak  for  N-H  and  O-H  protons  is  expected  to  occur  in  the  region  between  5  ~  6-8  ppm.  However, 
extensive  broadening  makes  peak  resolution  difficult  in  the  1 : 1  chloride  salt  case.  For  the  dication,  the 
peak  is  less  broadened  and  accurate  for  expected  integrations  of  each  proton  signal. 


600 


ppm(t1) 


9.0 


8.0 


7.0 


6.0 


5.0 


4.0 


3.0 


500 

400 

300 

200 

100 

0 

-100 


Figure  A4:  'H-NMR  plot  of  2-hydroxyethyl  hydrazinium  mono(tetrafluoroborate)  (top)  and  di(tetrafluoroborate)  (bottom) 
salts;  DMSO-</6  360  MHz. 


♦  There  is  an  unexplained  downfield  peak  at  5  ~  9.5  ppm,  and  the  peak  for  the  monocationic  protons 
integrates  to  less  than  expected  (2.75H(observed)  vs.  5H(expected)) 

♦  1:2  mixture  for  dicationic  salt  did  not  experience  the  above  problems,  and  the  obtained  sample,  however, 
the  proton  signal  for  -OH  and  ammonium/amine  protons  appears  to  be  somewhat  higher  in  integration 
than  expected  (~  9H  vs.  7H  expected).  Again,  this  might  be  attributed  to  the  difficulty  when  integrating 
the  broadened  signal. 
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Figure  A5:  'H-NMR  plot  of  2-hydroxyethyl  hydrazinium  mono(hexafluorophosphate)  (top)  and  di(hexafluorophosphate) 
(bottom)  salts;  DMSO-</6  360  MHz. 


♦  The  peaks  for  amine/ammonium  and  hydroxy  protons  show  a  reverse  of  expected  chemical  shift  for  the 
hexafluorophosphate  salts  compared  with  the  other  spectra  (8m0nocation  =  7.896  ppm  vs.  Station  =  6.953 
ppm). 

Part  B:  Thermal  Analysis  (TGA  / DSC  /  Thermal  Stability): 

♦  The  thermal  properties  for  the  synthesized  2-OHEt-HY  salts  are  summarized  in  the  table  below  (Table 
IB); 


Table  IB:  Thermal  Data  Collected  for  Mono-  and  Dicationic  2-OHEt  Hydrazinium  Salts 


SALT/RATIO 

TGA  (°C)t 

DSC  (°C)* 

Thermal  Stablility* 

2 -hydroxy  ethyl- 
hydrazine  (neutral) 

Tonset  =  181.98 

T  5%  dec  =  102.10 

Tdec=  195.94 

Tg  =  -72.1 1 

-  88.1  %  weight  loss 

2-OHEt  hydrazinium 
nitrate;  1:1 

Tonset  =  214.10 

T  5%  dec  =  193.63 

T  dee  =  259.67 

Tg  = -56.91 

-23.8  %  weight  loss 

2-OHEt  hydrazinium 
nitrate;  1:2 

Tonset  =91.51 

T  5%  dec  =  62.71 

T  dee  =  101.01 

Tg  =  -47.89 

-43.8%  weight  loss 

2-OHEt  hydrazinium 
chloride;  1:1 

TonSet  =  223.19 

T  5%  dec  =  171.15 

T  dec  =  25 1.33 

Tg  = -50.61 

-8.8  %  weight  loss 

2-OHEt  hydrazinium 
chloride;  1:2 

Tonset  =  240.33 

T  5%  dec  =  183.17 

Tg  = -49.37 

-7.8  %  weight  loss 

SALT/RATIO 

TGA  (°C)t 

DSC  (“C)* 

Thermal  Stablility* 

T  dec  =  282. 18 

2-OHEt  hydrazinium 
tetrafluoroborate;  1:1 

Tonset  =  198.49 

T  5%  dec  =  192.79 

T  dee  =  310.59 

Tg  =  -46.22 

-7.1  %  weight  loss 

2-OHEt  hydrazinium 
tetrafluoroborate;  1:2 

Tonset  =  160.66 

T  5%  dec  =  79.02 

T  dec  =  215.26 

Tg  = -55.49 

-38.8  %  weight  loss 

2-OHEt  hydrazinium 
hexafluorophosphate; 

1:1 

Tonset  =  106.56 

T  5%  dec  =  1  19-21 

T  dec  =  114.67 

Tg  = -32.75 

-40.2%  weight  loss 
(-29.2  %  weight  loss 
with  sample  dried  at  60 
°C  for  12  hours) 

2-OHEt  hydrazinium 
hexafluorophosphate; 

1:2 

Tonset  =  106.91 

T  5%  dec  =  1  19-21 

T  dec  =  114.67 

Tg  = -33.46 

-  41.5  %  weight  loss 

2-OHEt  hydrazinium 
bis(trifylamide);  1:1 

Tonset  =  162.57 

T  5%  dec  =  143.75 

T  dec  =  183.03 

Tg  = -48.75 

-17.5  %  weight  loss 

2-OHEt  hydrazinium 
bis(trifylamide);  1:2 

Tonset  =  158.12 

T  5%  dec  =  156.12 

T  dec  =  178.34 

Tg  =  -51.04 

-12.2  %  weight  loss 

1  All  TGA  measurements  were  made  on  dynamic  heating  regime  using  a  TGA  model  2950  TA  Instruments.  Instrument  under  air 
atmosphere  and  samples  limited  to  below  15  mg.  Samples  were  first  heated  to  75  °C  and  kept  at  this  temperature  for  30  min  to  allow 
all  remaining  water  to  evaporate.  Then  the  samples  were  heated  at  a  rate  of  5  °C  min'1  from  room  temperature  to  600  or  800  °C. 

‘  DSC  analysis  was  conducted  on  a  TA  Instruments  Modulated  DSC  cooled  with  a  liquid  nitrogen  cryostat,  and  data  was  collected  at 
atmospheric  pressure  using  5-20  mg  samples  in  sealed,  pin-holed  aluminum  pans  with  caps.  An  empty  pan  was  used  as  reference. 

*  Thermal  stability  tests  were  taken  for  samples  by  using  -10-20  mg  samples  with  the  TGA  and  setting  to  ramp  at  5  °C  min"1  to  75  °C 
and  allowing  isothermal  for  24  hours. 
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Figure  IB:  Comparative  TGA  plot  for  the  thermal  decomposition  of  monocation  (dotted  lines)  and  dicationic  (solid  lines)  salts 


of  2-hydroxyethyl  hydrazine  (black  line,  reference). 


Thermal  decomposition :  Thermograviometric  analysis  (TGA)  was  utilized  to  identify  the  temperatures  for 
decomposition  (Tdec)  and  temperature  for  onset  of  decomposition  (TonSet).  In  addition,  the  temperature  at  which 


5  %  of  the  material  has  thermally  decomposed  (Tso/0  dec)  is  included  as  a  more  accurate  assessment  of  the  real 
onset  for  thermal  decomposition  of  the  bulk  material.  Conditions  for  these  analyses  were  to  heat  5-15  mg  of 
sample  to  800  °C.  The  TGA  plots  (Figure  IB)  provide  a  comparative  view  of  how  the  series  of  synthesized  salts 
behaved  when  undergoing  thermal  decomposition. 

Thermal  decomposition  results  discussion: 

♦  Generally,  the  monocation  for  nitrate  seems  to  be  more  stable  than  the  dication  as  can  be  seen  when 
noting  the  earlier  onset  of  decomposition  for  the  dication. 

♦  The  ratio  of  mass  lost  in  the  first  and  second  decomposition  steps  of  the  dicationic  nitrate  might  suggest 
the  loss  of  (2)  molecules  of  HNO3. 

♦  In  general,  the  chloride  salts  (mono-  and  dicationic  forms)  demonstrated  very  good  thermal  stabilities, 
with  quite  a  late  onset  of  decomposition. 

♦  Additionally,  the  decomposition  of  the  chloride  salts  occurred  with  the  majority  of  weight  lost  in  the 
first  step.  A  minor  decomposition  step  can  be  observed  near  450  °C. 

♦  The  BF4'  salts  both  demonstrated  a  fairly  direct  decomposition  near  200  °C,  however. 

♦  The  explanation  for  differences  in  the  onset  of  decomposition  may  be  due  to  changes  in  the  composition 
of  the  salt  mixture  that  result  from  one  or  more  of  the  following  reason: 

o  Change  in  composition  as  a  result  of  dynamic  acid/base  equilibrium  where  only  one  of  the 
components  present  (acid  or  starting  material)  evaporates. 

♦  PF6’  salts  were  found  to  degrade  immediately  upon  the  start  of  TGA,  even  without  the  presence  of  an 
isothermal  event. 

♦  Both  of  the  NTf)"  salts  were  found  to  have  (3)  distinct  decomposition  events  involved  with  the  loss  of 
weight. 
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Figure  2B:  Comparative  TGA  plots  for  the  thermal  stability  of  obtained  monocationic  (dotted  lines)  vs.  dicationic  (solid  lines) 

salts  of  2-hydroxyethyl  hydrzazine  (black  line,  reference). 

Thermal  Stability.  The  relative  thermal  stabilities  of  each  salt  is  assessed  by  isothermic  heating  of  a  small 
sample  on  the  TGA  instrument  for  24  hours  @  75  °C.  The  %  weight  loss  of  the  sample  is  provided  above  in 
Table  Bl,  and  the  comparative  TGA  plots  for  thermal  stability  are  provided  in  Figure  2B.  Conditions  for  theses 
analyses  were  as  follows: 


o  Ramp  of  temperature  by  5  °C/min  to  75  °C. 
o  Isothermal  heating  for  24  hours. 


Thermal  Stability  Data  Discussion: 


♦  The  presence  of  an  initial  onset  that  varies  to  an  arbitrary  degree  might  suggest  the  possibility  of 
dynamically  changing  composition  due  to  the  loss  of  either  acid  or  starting  material  from  evaporation 
(as  mentioned  earlier). 

♦  Such  confounding  data  beckons  further  investigation  to  probe  into  the  cause  for  the  early  onset  of 
decomposition  in  the  isothermal  data  as  collected. 

DSC  Data  (based  on  data  summarized  in  Table  Bl): 

♦  Glass  transitions  are  generally  at  higher  temperatures  compared  with  the  parent  HEH. 

♦  All  salts  were  recorded  as  room  temperature  ionic  liquids  with  glass  transitions  <0  °C 

A.  Preliminary  Conclusions  and  Future  Directions: 

♦  Work  completed  includes: 

o  Synthesis  of  (10)  HEH-based  salts  on  1 .5  gram  scale  (Five  1 : 1  and  five  1 :2) 

■  Nitrates 

■  Chlorides 

■  Tetrafluoroborates 

■  Hexafluorophosphates 

■  Bistrifylamides 

o  Synthesis  of  2-OHEt-HY  [NO3]  1 : 1  salt  on  5  gram  scale 
o  Characterization: 

■  'H-NMR 

■  TGA 

■  DSC 

■  Thermal  Stability  Testing 

♦  The  prototype  of  the  1 : 1  and  1 :2  nitrate  salts  indicate  higher  thermal  stability  of  the  salts  over  the  neutral 
HEH  compound;  TGA  analysis  indicates  a  more  stable  monocationic  salt  when  compared  with  the 
dicationic  dinitrate. 

♦  Potential  anion  effects  exist  as  well  as  speciation  on  phase  transitions,  as  noted  in  the  DSC  data  requires 
further  attention. 
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